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We measure the branching fraction and the CP-violating asymmetry of B° —> K° n° decays with 
227 million T(4S) —> BB events collected with the BABAR detector at the PEP-II asymmetric-energy 
e + e _ collider at SLAC. We obtain a branching fraction B(B° —> K°n 0 ) — (11.4 dt 0.9 ± 0.6) x 10~ 5 
and CP-violating asymmetry parameters C K o n o = 0.06 ± 0.18 ± 0.03 and S K o,o = 0.351^33 ±0.04, 
where the first error is statistical and the second systematic. 

PACS numbers: 13.25.Hw, 13.25.-k, 14.40.Nd 


CP violation effects in decays of B mesons that are 
dominated by penguin b —> sqq (q = u, d, s ) tran¬ 
sitions are potentially sensitive to contributions from 
physics beyond the Standard Model (SM) [1], The B- 
factory experiments have explored time-dependent CP- 
violating (CPV) asymmetries in several such decays [2], 
namely B° -> <j>K° s [3, 4], B° -> rfK° [3, 5], B° -► 
K+K~K° S [3, 6 ], B° -> f 0 K° s [7] and B° -» K° s ir° [ 8 ], 
Within the SM these asymmetries are expected to be 
consistent with the measurement of sin2/3 in charmonium 
modes originating from the tree-level b —> ccs transition. 
These comparisons must take into account possible devi¬ 
ations for each mode, within the SM, due to contributions 
of other diagrams with different phases and rescattering 
effects. At this point none of the modes above shows a 
significant deviation from the SM expectation [9]. A ma¬ 
jor goal of the P-factory experiments is to reduce the ex¬ 
perimental uncertainties of these measurements in order 
to improve the sensitivity to beyond-the-Standard-Model 
effects. 

In this letter we present improved measurements of 
the CPV asymmetry in the decay B° —> using 

data collected with the BABAR detector at the PEP-II 
asynnnetric-energy e + e _ collider, amounting to 226.6 ± 
2.5 million T(4S ) —» BB decays. In the SM this de¬ 
cay is dominated by a top-quark-mediated b —> sdd pen¬ 
guin amplitude. If other contributions, such as the CKM 
suppressed b — > suu tree amplitude, are ignored, the 
CPV asymmetry is governed by sin2/3 [10], where f3 = 
arg [—V c dV* b /VtdV t * b ] and V is the Cabibbo-Kobayashi- 
Maskawa (CKM) quark mixing matrix [11]. The bound 
on the deviation from sin2/I due to SM contributions with 
a different weak phase is about 0.2 from SU(3) flavor sym¬ 
metry [12] and about 0.1 in model-dependent QCD cal¬ 
culations [13]. We also present an update of our measure¬ 
ment of the branching fraction of B° —» K°tt° [14], which, 
when combined with measurements of other B —> Kn 
branching fractions, can be used to extract the CKM an¬ 
gle 7 = arg[-V ud V* b /V cd V* b \ [15]. 

The BABAR detector, fully described in [16], pro¬ 
vides charged particle tracking through a combination 
of a five-layer double-sided silicon micro-strip detector 
(SVT) and a 40-layer drift chamber (DCH), both oper¬ 
ating in a 1.5 T magnetic field. Charged-kaon and -pion 
identification is achieved through measurements of spe- 
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cific energy-loss ( dE/dx) in the tracking system and of 
the Cherenkov angle (9 C ) in a detector of internally re¬ 
flected Cherenkov light (DIRC). A CsI(Tl) electromag¬ 
netic calorimeter (EMC) provides photon detection and 
electron identification. Finally, the instrumented flux re¬ 
turn (IFR) of the magnet allows discrimination of muons 
from pions. For event simulation we use the Monte Carlo 
event generator EVTGEN [17] and GEANT4 [18]. 

At the T(4S) resonance time-dependent CPV asymme¬ 
tries are extracted from the distribution of the difference 
of the proper decay times, A t = top — ftag> where tcp 
refers to the decay time of the signal B ( Bcp ) and f ta g 
to that of the other B (B ta , g ). The At distribution for 
Bcp f follows 

— |At|/r 

V ± (At) = ( 1 ) 

[1 ± Sf sin (AmrfAt) Cf cos (Am<jAt) ] , 

where the upper (lower) sign corresponds to B ta , g decay¬ 
ing as B° (B°), t is the B° lifetime and A rrid is the 
mixing frequency. The coefficients Cf and Sf can be ex¬ 
pressed in terms of the B°-B° mixing amplitude and the 
decay amplitudes for B° —» / and B° —> / [19]. For de¬ 
cays to a CP eigenstate, like K® 7 r°, Cf vanishes unless 
there is direct CP violation. If B° —» K®n° proceeds 
purely through a top-quark penguin, C K o v o = 0 and 
S K o n o = sin(2/3 + 2 f3 s ), where /3 S = a,rg[-V ts V* b /V cs V* b ] 
is small. 

We search for B° —» A"° 7 r° decays in BB candidate 
events selected using charged-particle multiplicity and 
event topology [20]. We reconstruct K° —> candi¬ 

dates from pairs of oppositely charged tracks. The two- 
track combinations must form a vertex with a y 2 con¬ 
sistency greater than 0.001 and a 7 r + 7 r - invariant mass 
within 11.2 MeV/c 2 of the nominal K® mass [21], We 
form 7 r° —» 77 candidates from pairs of photon candi¬ 
dates in the EMC, each of which is isolated from any 
charged tracks, carries a minimum energy of 50 MeV, 
and has the expected lateral shower shape. Candidates 
for B° —> A'° 7 t° are formed from K®n° combinations and 
constrained to originate from the e + e _ interaction point 
using a geometric fit. We require that the consistency 
of the x 2 °f the fit: which has one degree of freedom, be 
greater than 0.001. We extract the AT° decay length L K 0 
and the 7r° —> 77 invariant mass from this fit and require 
110 < m 77 < 160 MeV/c 2 and L K 0 greater than 5 times 
its uncertainty. 

For each B candidate we compute two kinematic vari¬ 
ables, namely the invariant mass mp and the missing 
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mass TO m i ss = \J(q e + e - — <}b) 2 , where q e + e - is the four- 
momentum of the initial e + e~ system and qs is the 
four-momentum of the B° —> K®tt q candidate after a 
mass constraint on the B° is applied. By construction 
the linear correlation coefficient between ?n m i ss and mp 
vanishes. Compared to the kinematic variables A E = 

e *b ~ \\fs an d ?ties = \J\s - p*p (where s = g 2 +e _ 
and the asterisk denotes the e + e^ rest frame), which 
were used in our previous analysis of this mode [8], the 
present combination of variables leads to a smaller cor¬ 
relation and a better background suppression for modes 
containing a high-momentum ir° or photon. From sim¬ 
ulation studies we determine the signal resolution for 
m_e to be about 40 MeV/c 2 . The distribution exhibits 
a low-side tail from energy leakage out of the EMC. 
The signal resolution for m m ; ss , about 5 MeV/c 2 , is dom¬ 
inated by the beam-energy spread. We select candi¬ 
dates with rriB within 150 MeV/c 2 of the nominal B° 
mass [21] and with 5.11 < m m i ss < 5.31 GeV/c 2 . The 
region m m ; ss < 5.2 GeV/c 2 is devoid of signal and used 
for background characterization. 

To suppress background from continuum e + e~ —* qq 
(■q = u, d , s, c) events, we exploit differences in both pro¬ 
duction and decay properties. We require | cosd^| < 0.9, 
where 0* B is the angle between the B-candidate momen¬ 
tum and the e~ momentum in the e + e _ rest frame. For 
true B mesons the distribution of cos 0* B is proportional 
to 1 — cos 2 9* b , whereas for continuum events it is nearly 
flat. To exploit the jet-like topology of continuum events, 
we calculate the ratio L 2 /B 0 of two Legendre moments 
defined as Lj = JV |p*|| cosd*p, where p* is the momen¬ 
tum of particle i in the e + e~ rest frame, 9* is the angle 
between p* and the thrust axis of the B candidate and 
the sum runs over all reconstructed particles except for 
the B-candidate daughters. We require L 2 /Bo < 0.55, 
which suppresses the background by more than a factor 
3 at the cost of approximately 10% loss in signal effi¬ 
ciency. After all selections are applied the average can¬ 
didate multiplicity in events with at least one candidate 
is approximately 1.007. When there are multiple candi¬ 
dates, we select the candidate with a reconstructed 7r° 
mass closest to the expected value. 

For each B° —> K°7r° candidate we examine the re¬ 
maining tracks in the event to determine the decay ver¬ 
tex position and the flavor of B tag . Using a neural 
network based on kinematic and particle identification 
information [22] each event is assigned to one of five 
mutually exclusive tagging categories, designed to com¬ 
bine flavor tags with similar performance and At res¬ 
olution. We parameterize the performance of this al¬ 
gorithm in a data sample (Bfl av ) of fully reconstructed 
B° —» D < -*' > ~n + / p + /af decays. The average effective 
tagging efficiency obtained from this sample is Q = 
E c e|(l - 2u; c ) 2 ' = 0.299 ± 0.005, where e c s and w c are 
the efficiencies and rnistag probabilities, respectively, for 
events tagged in category c = 1, 2, • • • 5. For the back¬ 
ground, the fraction of events ( e B ) and the asymmetry 


in the rate of B° versus B° tags in each tagging category 
are extracted from a fit to the data. 

The proper-time difference is extracted from the sep¬ 
aration of the Bcp and B tag decay vertices. The B tag 
vertex is reconstructed inclusively from the remaining 
charged particles in the event [20]. To reconstruct the 
Bqp vertex from the single K° trajectory we exploit the 
knowledge of the average interaction point (IP), which is 
determined on a run-by-run basis from the spatial dis¬ 
tribution of vertices from two-track events. We com¬ 
pute At and its uncertainty from a geometric fit to the 
T(4S) —> B°B° system that takes this IP constraint into 
account. We further improve the sensitivity to At by con¬ 
straining the sum of the two B decay times (tcp+t t ag ) to 
be equal to 2 r B 0 with an uncertainty v /2 rpo , which ef¬ 
fectively constrains the two vertices to be near the T(4S) 
line of flight. We have verified in a Monte Carlo simu¬ 
lation that this procedure provides an unbiased estimate 
of At. 

The per-event estimate of the uncertainty on At re¬ 
flects the strong dependence of the At resolution on the 
Kg flight direction and on the number of SVT layers tra¬ 
versed by the K° decay daughters. In about 60 % of the 
events both pion tracks are reconstructed from at least 
4 SVT hits, leading to sufficient resolution for the time- 
dependent measurement. The average At resolution in 
these events is about 1.0 ps. For events which fail this 
criterion or for which cr(Af) > 2.5 ps or At > 20 ps, the 
At information is not used. However, since Cf can also 
be extracted from flavor tagging information alone, these 
events still contribute to the measurement of Cf. 

We extract the signal yield, Sf and Cf from an un¬ 
binned maximum-likelihood fit to ms, m m j ss , L 2 /L 0 , 
cos 9 b , At and the flavor tag variables. By exploiting 
sideband regions in data for the background and simu¬ 
lated events for the signal, we have verified that with 
the selection presented above these observables are suffi¬ 
ciently independent that we can construct the likelihood 
from the product of one-dimensional probability density 
functions (PDFs). The PDFs for signal events are param¬ 
eterized from simulated events or from the Bfl av sample. 
For background PDFs we select a functional form that 
describes the data in the sideband regions of the other 
observables, in which backgrounds dominate. We include 
these regions in the fitted sample and simultaneously ex¬ 
tract the parameters of the background PDFs along with 
the signal yield and CPV asymmetries. 

We obtain the PDF for the A t of signal events from the 
convolution of Eq. 1 with a resolution function lZ(5t. = 
At — Af t rue 5 crAt)- The resolution function is parameter¬ 
ized as the sum of two Gaussians with a width propor¬ 
tional to the reconstructed a At, and a third Gaussian 
with a fixed width of 8 ps [20]. The first two Gaussians 
have a non-zero mean, proportional to o At, to account 
for the small bias in At from charm decays on the B tag 
side. We have verified in simulation that the parame¬ 
ters of TZ(8t, UAt) for B° —> A'g 7 r° events are similar to 
those obtained from the Bfl av sample, even though the 
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FIG. 1: Signal and background (inset) distributions, obtained 
with the weighting technique described in the text, for ms 
(left) and m m i ss (right). The curves represent the PDFs used 
in the fit and are normalized to the fitted yield. 

distributions of a At differ considerably. We therefore ex¬ 
tract these parameters from a fit to the Bfl av sample. We 
assume that the background consists of prompt decays 
only and find that the At distribution is well described 
by a resolution function with the same functional form as 
used for signal events. The parameters of the background 
function are determined in the fit. 

To extract the yield and the CPV asymmetries we 
maximize the logarithm of the extended likelihood 

C(S f ,Cf,Ns,N B ,fs,fB,d) = e~ (Ns+NB) x 

[Nsfs£s'Ps(xi,yi-,Sf,Cf) + N B fBe C B'PB(xi,yi;d)] x 
iei 

n i N s( 1 - fsVsV's (ft; Cf) +Nb( 1 - a)] , 

ieii 

where I (II) is the subset of events with (without) 
At information. The probabilities Vs (V' s ) and Vb 
(V' b ) are products of PDFs for signal ( S) and back¬ 
ground ( B) hypotheses evaluated for the measurements 
Xi = {m B , m miss ,L 2 /io, cos 0* B , tag, tagging category} 
and j/i = {At, a At}- Along with the signal yield N$ 
and the coefficients Sf and Cf , the fit extracts the back¬ 
ground yields N B , the fractions of events with At infor¬ 
mation fs and f B , and the remaining parameters, col¬ 
lectively denoted by a. These include all parameters 
of background PDFs and some parameters of the signal 
PDFs, such as the mean values of m B and m m i ss . 

Fitting the data sample of 9726 B° —> K®tt° can¬ 
didates, we find Ns = 300 ± 23 signal decays with 
S K o n o = 0-35 io .33 (stat) ± 0.04 (syst) and C K o^o = 
0.06 ± 0.18 (stat) ± 0.03 (syst). The number of signal 
decays with At information is f s N s = 186 ± 18. The 
total detection efficiency for B° —> K®7 r° decays with 
K° —> 7 r + 7 T - and 7 r° —► 77 is (34.1 ± 1.8) %. With the K® 
and 7T° branching fractions taken from [21] and assum¬ 
ing equal production of charged and neutral B mesons 
at the T(4S) resonance, we obtain a branching fraction 
B(B° -> K°ir°) = (11.4±0.9(stat)±0.6(syst)) x 10 " 6 . 
The evaluation of the systematic uncertainties is de¬ 
scribed below. 

Figure 1 shows the background-subtracted distribu- 



At(ps) 

FIG. 2: Signal distribution for At, obtained with the weight¬ 
ing technique described in the text, with Stag tagged as B° 
(top) or B° (center), and the asymmetry A K o n o(At) (bot¬ 
tom). The curves represent the PDFs for signal decays in the 
likelihood fit. 


tions oi m B and m m i ss for all B° —> K°ir° candidates in 
the fit. The background subtraction is performed with an 
event weighting technique [23]. Events contribute accord¬ 
ing to a weight constructed from the covariance matrix 
for the yields (Ns and N B ) and the probability Vs and 
V B for the event, computed without the use of the vari¬ 
able that is being displayed. The curves represent the sig¬ 
nal PDFs used in the fit. The insets show the correspond¬ 
ing signal-subtracted distributions with the background 
PDFs. Figure 2 shows the background-subtracted distri¬ 
butions of At for B°- and l?°-tagged events, and of the 
asymmetry A K o n o(At) = [N B 0 - N Bo ] / [N B o + N Bo ] as 
a function of At. 

The extraction of At with the IP-constrained fit has 
been extensively tested on large samples of simulated 
B° — > K®tt 0 decays with different values of C and S. 
We have also exploited a control sample of approximately 
1900 observed B° —> J/ipK^ decays with J/ij) —> YY 
and J/il> —> e + e _ , using the procedure described in [ 8 ]. 
Based on these studies we assign a systematic uncertainty 
of 0.023 on S and 0.014 on C due to the At reconstruc¬ 
tion and the choice of the resolution function. As a cross¬ 
check we measure the B° lifetime in B° —> decays 

in data and find that it agrees with the world average. 
We evaluate the effect of a possible misalignment of the 
SVT by introducing misalignments in the simulation and 
assign a systematic uncertainty of 0.020 on S and 0.007 
on C. We also consider large variations of the position 
and size of the interaction region, which we find to have 
negligible impact. We include a systematic uncertainty 
of 0.012 on S and 0.018 on C to account for imperfect 
knowledge of the PDFs used in the fit. Using simulated 


































7 


events we estimate a contribution of 2.3± 1.7 events from 
other B decays for which we assign a systematic uncer¬ 
tainty of 0.019 on S and 0.015 on C. 

The detection efficiency for signal events is obtained 
from a Monte Carlo simulation. The efficiency of the K ° 
selection is calibrated with a large sample of inclusive 
K® —> 7 r + 7 r~ decays. The n 0 —> 77 efficiency is cali¬ 
brated with e + e _ —> t + t~ events with r~ —* p~i' T . The 
systematic uncertainty associated with the efficiency is 
2.8% for K° s and 3.0% for 7 r°. We assign additional sys¬ 
tematic uncertainties of 1.2 % for the L 2 /B 0 cut, 2.0 % for 
the selection on ms and a total of 2 . 0 % for uncertain¬ 
ties in the signal PDFs. Finally, we include a systematic 
uncertainty of 1.4% to account for unknown contribu¬ 
tions from other B decays and a systematic uncertainty 
of 0 . 6 % due to the uncertainty in the total number of 
T(4S) —> BB decays. 

In summary, we have reported improved measurements 
of the branching fraction and CP-violating asymmetry 
for the decay B° —> K° 7 r°. The measured values of S K o w o 


and C K Ojo are consistent with the Standard Model pre¬ 
dictions. The measured branching fraction is consis¬ 
tent with measurements from other experiments [24]. 
These results supersede our previous measurements of 
the branching fraction [14] and CPV asymmetries [ 8 ], 
which were based on a subset of the data presented here. 
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